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Abstract 

The temperature and polarization anisotropies of the Cosmic Microwave Background are 
analyzed under the hypothesis that the same inflationary seed accounting for protogalac- 
tic magnetism also affects the Einstein-Boltzmann hierarchy whose initial conditions are 
assigned for typical correlation scales larger than the Hubble radius after matter-radiation 
equality but before decoupling. Since the primordial gauge spectrum depends on a combi- 
nation of pivotal parameters of the concordance model, the angular power spectra of the 
temperature and of the polarization are computed, for the flrst time, in the presence of 
a putative large-scale magnetic fleld of inflationary origin and without supplementary hy- 
potheses. 
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1 Bootstrapping large-scale magnetism 



More than sixty years ago, an intense debate on the origin of galactic cosmic rays [1, 2J 
and the first ambiguous evidences of starhght polarization [3l H] led Fermi to propose the 
existence of a galactic magnetic field with approximate strength in the /iG rang^. In the 
same period it was correctly argued that cosmic rays above 10^^ eV were galactic (as opposed 
to solar) and the high degree of isotropy in the arrival directions could be explained by the 
presence of an irregular magnetic field of /xG strength able to scramble the trajectories of 
charged species (see, for instance, J3] for a terse account of this problem). 

Today we know, with a plausible degree of confidence, that the galactic field contains a 
large-scale regular component and a small-scale turbulent one, both having a local strength 
of a few fiG (see, for instance, [6]). While the turbulent component dominates in strength 
by a factor of a few, the regular component (with approximate correlation scale between 10 
to 30 kpc) imprints dominant drift motions as soon as the Larmor radius of cosmic rays is 
larger than the maximal scale of the (kinetic or magnetic) turbulences which is (9(100 pc). 
Clusters (i.e. gravitationally bound systems of galaxies) have been shown to possess a 
large-scale magnetic field, in the /iG range (more specifically between 10 and 100 nG) and 
with correlation scale larger than in the galactic case suggesting a quasi-flat magnetic power 
spectrum [7]. Superclusters (i.e. loosely bound systems of clusters) have been also claimed to 
have magnetic fields [,8j at the fiG level even if, in this case, unresolved ambiguities persist on 
the way the magnetic field strengths are inferred from the Faraday rotation measurements. 

In a seemingly different perspective we are now facing a problem conceptually close 
to the one of the early fifties of the past century. The spectrum of the highest energy 
cosmic rays shows no significant deviation from isotropy [9j as it can be argued from the 
distribution of arrival directions of cosmic rays detected above 10^^ eV at the Pierre Auger 
Observatory. Purported correlations between the arrival directions of cosmic rays with energy 
above 6 x 10^^ eV and the positions of active galactic nuclei within 75 Mpc [10] are then 
statistically insignificant pj. At smaller energies it has been convincingly demonstrated [TU] 
that overdensities on windows of 5 deg radius (and for energies 10^^'^eV < E < 10^^'^eV) are 
compatible with an isotropic distribution. We can then conclude that in the highest energy 
domain (i.e. energies larger than about 60 EeV), cosmic rays are not appreciably deflected 
and basically isotropic. As in the the early fifties the isotropy of the cosmic ray spectrum 
was a simple (thought indirect) evidence of galactic magnetism, it is tempting to speculate 
that the isotropy of the highest energy cosmic ray tells a similar story on irregular magnetic 
fields with correlation scales (9(/ewMpc) and within a cocoon of 70 Mpc. 

The current pieces of evidence confirm the early suggestions of Fermi arguing that large- 
scale magnetization should somehow be the remnant of the initial conditions of the gravi- 
tational collapse of the protogalaxy occurring over typical length-scales (9(Mpc). Harrison 



^We shall employ hereunder the conventional prefixes of International System of units i.e. /xG 
10"^ Gauss, nG = 10^^ Gauss and so on and so forth. 
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[12] and, with slightly different accents, Zeldovich [13] were among the first ones to pro- 
pose that the source of magnetism over the largest length-scales could be searched in the 
pre-decoupling plasma and, since then, various mechanisms have been tailored for the gener- 
ation of large-scale magnetic fields (see, for instance, [H]). This class of problems has been 
dubbed some time ago magnetogenesis [IH], i.e. the generation of large-scale magnetic fields 
during the early stages of the evolution of the plasma. Since the magnetic fields must not 
jeopardize the spatial isotropy of the background geometry, their Fourier modes must be 
stochastically distributed and characterized by an appropriate power spectrum: 

r) B,{k, r)) = ^ Ps{k, r) P;,{k) 5^'\q + k), (1.1) 

where, following the standard conventions (see e.g. [E]), Pis{k,T) is the physical power 
spectrum, Pij{k) = {6ij — kikj) and ki = ki/\k\. As it can be easily verified from the 
definition of the Fourier transform, P^ik^T) has dimensions of an energy density and its 
square root has, therefore, the dimensions of a field intensity. Since magnetic fields exist 
over increasing length-scales with almost the same intensity, the parametrization of Eq. 
f ll.ip implies that P^ik^r) can be considered, in the first approximation, as nearly scale- 
invariant. Incidentally, on a technical ground, the parametrization of Eq. (11. ip follows the 
same conventions employed when the power spectra of curvature perturbations are assigned 
at the standard pivot scale fcp = 0.002 Mpc~^; denoting with 7^(/c,r) ~ Tl^{k) the constant 
mode of curvature perturbations on comoving orthogonal hypersurfaces prior to equality, its 
two point function in Fourier space is 

(7^(g-)7^(^)) = ^PT^(k)6^'\q + k), P^^^k) = An[^y~\ (1.2) 

where, using the WMAP 9yr data alone [161 113 HH] in the light of the concordance scenario 
An = (2.41 ± 0.10) X 10"^ and = 0.972 ± 0.013. As in the magnetic case, the exact 
scale- invariant limit is realized when — )■ 1. 

Magnetic fields can be produced during a phase of decelerated expansion inside the 
particle horizon which always grows much faster than the correlation scale of the field. It is 
then unlikely to obtain, in this case, correlation scales (9(Mpc) at the onset of protogalactic 
collapse: this is one of the main drawbacks of what is often called, somehow improperly, 
causal magnetogenesis [IH] . Conversely during inflation the Hubble radius is almost constant 
0{H~^) and it roughly coincides with the event horizon. The correlation scale of a putative 
magnetic field amplified during inflation evolves much more rapidly than the Hubble radius 
itself so that the magnetic power spectrum at the end of inflation has an approximate 
amplitude 0{H'^). In the context of inflationary models the correlation length-scales of the 
produced magnetic fields range between 0{k~^) and 0{k^^) where k^ (see Eq. (11. 2p ) is 
the comoving scale at which the power spectrum of curvature perturbations is assigned in 
standard analyses and kj^ ~ (9(Mpc~^) corresponds to the wavenumber of the gravitational 
collapse of the protogalaxy. 
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To make any statement concerning the early evolution of the plasma, the post-inflationary 
dynamics must be specifled within a reasonable degree of accuracy. Direct Cosmic Microwave 
Background observations [HI [T71 HI] (CMB in what follows), light curves of type-la super- 
novae [201 EI] and large-scale galaxy surveys [22] are compatible with a wide class of scenarios 
arranged around the ACDM paradigm where A stands for the dark energy component and 
CDM for the cold dark matter component. The ACDM paradigm consists, nominally, of 6 
independent parameters which are extracted from the observational data. However, deflnite 
statements demand also more stringent assumptions on the origin of the adiabatic mode or 
on the thermal history of the plasma at least up to energy scales as large as few TeV. A 
conventional completion of the ACDM paradigm consists of a phase of slow-roll inflation 
with standard thermal history where the post-inflationary evolution is suddenly dominated 
by radiation down to the scale of matter-radiation equality. In the latter framework the 
CMB scale left the inflationary Hubble radius around 65 efolds prior to the end of inflation; 
conversely the scale of the protogalactic collapse left the inflationary Hubble radius around 
9 efolds after the CMB scale. 

The conventional completion of the ACDM paradigm implies that the typical wavenumber 
of the gravitational collapse of the protogalaxy in units of the inflationary Hubble rate is: 

^ = 3.22 X 10-- (J^) (^Y'^' ( y\ (1.3) 

aH VMpc-V VO.Ol/ V2.41 x lO'V ' ^ ' 

where, besides An, e = —H/H^ deflnes the standard slow-roll parameter. The quasi flat 
spectrum of inflationary origin at the present time and for a typical comoving scale 0{ki) 
can then be fully specifled in terms of the standard thermal history and it is given bjifl 

Ps(k,, r„) = 10-- r^^^l m {^) nG^ (1.4) 



,2.41 X 10-V VO.Ol/ V4.15 x 10-V V0.22 

where it has been assumed Psik, Tend) = rH'^ ~ 0(i7^) at the end of inflatioifl. The 
magnetic fleld itself can by estimated by taking the square root of the Eq. (11. 4p giving a 
fleld of the order of 0.06 nC Equation (ll.4p is compatible with a galactic magnetic fleld 
of the order of the /xG since during the process of collapse (prior to the onset of galactic 
rotation) the magnetic flux is frozen into the plasma element thanks to the large value of the 
conductivity. The mean matter density increases, during collapse, from its critical value (i.e. 
Per = 1-05 X 10^^ /igGeV/cm^) to a flnal value pf value which is 5 to 6 orders of magnitude 
larger than pc. The magnetic power spectrum after collapse will then be larger than (ll.4p 
by a factor which is roughly (pf / p, 



■^Note that /iq^r corresponds to the present critical fraction of radiation given as the sum of the energy 
density of the photons and of the neutrinos which are strictly massless in the ACDM scenario. 

*In specific models it turns out that the power spectrum, in the quasi-flat case, is given by r H'^ where 
r is a numerical factor ranging between 0.1 and 0.01. The reference value of r is r = 9/(47r^) ~ 0.22 and it 
corresponds to the one obtainable in the context of explicit models mentioned later on. 



4 



As a consequence of the standard thermal history, the amphtude of the magnetic power 
spectrum depends only on e and Atz and the inflationary Hubble rate in Planck units is 
given hj H = a/tt e Atz Mp . But since An and h^Q-R are fixed within the vanilla ACDM 
scenario, the amplitude of the seed spectrum effectively depends only on e. The purpose of 
this paper is to determine the initial conditions of the Einstein-Boltzmann hierarchy in terms 
of the early initial conditions derivable in the framework of inflationary magnetogenesis. It 
has been pointed out (see {23| and references therein) that a key role in the analysis of 
large-scale magnetism must rely on an accurate determination of the distortions induced 
on the temperature and polarization anisotropic^ and various results have been derived 
both analytically and numerically but never in terms of a primordial inflationary seed. The 
aim here is to bridge this gap by treating large-scale fields not generically present prior 
to decoupling but rather coming from inflationary magnetogenesis and to derive the initial 
conditions of the Einstein-Boltzmann hierarchy. The TT, EE and TE angular power spectra 
will then be directly computed in terms of the appropriate initial conditions. 

To develop a theory of the magnetized initial conditions without running into well known 
troubles, an efficient tool is the synchronous gauge which is particularly suitable for the 
perturbative description of the anisotropic stresses. The synchronous results will be explicitly 
connected to the set of gauge-invariant variables derived by Lukash [21] in the context of the 
Lifshitz formalism [22] (see also [5^1 EZ])- During the post-inflationary epoch, the evolution 
of the whole irrotational system is reduced to a single normal mode which is invariant under 
infinitesimal coordinate transformations (as required in the context of the Bardeen formalism 
[28] ) and whose source terms depend on the magnetic inhomogeneities. In the absence of 
magnetic fields, earlier analyses [291 EO] followed the same logic of [21] but in the case of scalar 
field matter. We shall use that the normal modes identified in [2H 123 EO] are all related to 
the (rescaled) curvature perturbations on comoving orthogonal hypersurfaces [511 E2] • 

The layout of the present analysis is as follows. In section[2]the fluctuations of inflationary 
magnetogenesis are analyzed in the synchronous coordinate system. In section |3] the initial 
conditions of the Einstein-Boltzmann hierarchy are bootstrapped out of the values provided 
by inflationary magnetogensis by using the explicit evolution of curvature perturbations for 
typical scales larger than the Hubble radius at the corresponding epoch. The whole system of 
perturbations is then integrated across the radiation-matter transition. In section H] the TT, 
EE and TE correlations are explicitly computed in the case of a magnetized adiabatic mode 
of inflationary origin. The B-mode autocorrelations arising in this approach are also briefly 
examined. Section O contains the concluding remarks. In the appendix some of technical 
results have been collected for the interested reader to avoid excessive digressions. 

^In the ACDM paradigm the temperature and polarization anisotropics are given in terms of the temper- 
ature autocorrelations (the TT power spectrum), the polarization autocorrelations (the EE power spectrum) 
and the temperature-polarization cross-correlations (the TE power spectrum). The other correlations are 
vanishing even if, in the presence of a stochastic magnetic field, a B-mode power spectrum is naturally 
induced by Faraday rotation. 
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2 Inflationary evolution 



During a phase of slow-roll 




magnetic fields are amplified while the electric fields 



are suppressed thanks to the coupling to a spectator field [331 EH [35]. The electric and 
magnetic power spectra contribute to the fluctuations of the various components of the 
energy-momentum tensor and, ultimately, to the curvature perturbations in the same way 
as a putative magnetic field affects curvature perturbations across the matter radiation 
transition (see [19] and references therein). The form of the electric and magnetic power 
spectra will now be discussed not with the aim of endorsing a particular magnetogenesis 
scenario but with the purpose of drawing general lessons on the parametrization of the 
magnetized power spectra of curvature perturbations at the end of inflation. 

2.1 Parametrization of the power spectra 

Consider an explicit coupling of the Abelian kinetic term to a spectator field a: 



where (p denotes the infiaton field, the gauge field, Y^y the gauge field strength; rep- 
resents the potential contribution of an Ohmic current. The evolution of the gauge fields in 
this class of models has been analyzed in different situations (see, e.g. [331 El] and references 
therein). If A (a) is a monotonically increasing function of the conformal time coordinate 
r, magnetic fields are amplified while the electric fields are suppressed. This kind of mono- 
tonic behavior is quite natural in conventional inflationary models where the curvature scale 
increases as we approach the protoinfiationary stage which is customarily modeled with a 
phase of decelerated expansion preceding the ordinary slow-roll phase. During the protoinfia- 
tionary phase conformal invariance can still be unbroken and, if this happens, initial Ohmic 
currents (remnants of the protoinfiationary dynamics) are not dissipated [3S]- Conducting 
initial conditions may then dominate against the quantum initial conditions and the increase 
of A corresponds, in this case, to an increase of the Debye shielding length and to a further 
suppression of the electric fields [35]. For an exponential coupling we shall have, for instance. 
In A (cr) = 7o-cr/M where M is a typical scale characterizing the evolution of the spectator 
field. The equations of motion of a can be solved at the level of the background and more 
specific situations can be found, for instance, in [33] (see also [31] )■ 

^Consistently with the absence of spatial curvature of the ACDM paradigm, the background geometry is 
assumed to be conformally flat throughout this investigation; the line element is given by ds^ = a? (r) [o?t^ — 
dx^]. The prime denotes a derivation with respect to the conformal time coordinate r while the overdot 
indicates a derivation with respect to the cosmic time coordinate t. The conformal time coordinate r is 
related to the cosmic time as a{T)dT = dt. 



s 




(2.1) 
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In a conformally flat background geometry characterized by a scale factor a(r) in the 
conformal time parametrization, the power spectra of the comoving electric and magnetic 
fields ar^: 

P^{k,T) = ^\MT)\', f,{r) = ^V^Hi'\-kT), (2.2) 




PE{k,T) = ^\gkir)\\ g,{r) = ^ -V-kr Hl'M-kr), (2.3) 

where the mode functions /^(t) and (7a,. (r) are solutions of the corresponding equation^: 

f'k = 9k + J^fk, 9'k = -k'^fk- J^9k - 47ro-c5ffc- (2.4) 

In Eq. (12.41) the prime denote a derivation with respect to the conformal time coordinate 
and T = / VX- To keep the discussion general without selecting a specific model, the 
growth rate J-" is parametrized as J-" = fH with T-L = aH; in the latter case the Bessel 
index appearing in Eqs. ( 12. 3p and ( 12. 2p isz/ = l/2 + / + /e and e is the slow-roll parameter 
entering the expression of u since the relation between the Hubble rate and the conformal 
time coordinate demands the following well known condition (1 — e)aH = — r. In the case 
of Eq. (12.11) / = e7o-/(l — e) but the parametrization in terms of the growth rate has the 
advantage of being sufficiently general to incorporate at once different dynamical situations. 
The quasi-flat magnetic spectrum corresponds to the case / ~ 2; direct analyses of these 
scenarios show that, in this class of models, the departure from scale-invariance cannot be too 
large (i.e. / < 2.2) [36] if the adiabatic mode is to be the dominant source of inhomogeneities 
across matter-radiation decoupling. After the end of inflation J-" — )■ and the relation 
between the physical and the comoving power spectra is given by P-Q{k,r) = a^{T)Pts{k,T) 
and PE(fc,r) = a'^{T)P£{k,T). 

As discussed in appendix lAl the fields B{x,t) and E{x,t) correspond to the canonical 
normal modes of the system diagonalizing the Hamiltonian density and having well defined 
properties of transformation under the duality symmetry [351 EH] • In terms of the power 
spectra of the comoving magnetic and electric fields the correlation functions are given by 

(S,(g,r)5,(p,r)) = ^ P^{q,r) P,^{q) 6^'\q + p), (2.5) 

{E,{q,r)Ej{p,T)) = ^ P^{q,r) P,,{q) 6^'\q + p). (2.6) 

In full analogy with the power spectra of curvature perturbations introduced in section [H 



^Note that H,y{z) are standard Hankel functions. Furthermore \M\'^ = 7r/2. 

®The presence of the conductivity CTc (not to be confused with the spectator field) plays a specific role 
during reheating and in the case of conducting initial conditions. More details on these issues can be found in 
[35j . Here we shall focus, for illustration, on the conventional case of quantum mechanical initial conditions, 
as it can be argued from the boundary conditions imposed on the mode functions. 
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the magnetic power spectra can be expressed as 



PB{k,T) = AB{K,T)[—] , 



9jr 

47]-2 



kLJ 
[aHj 



nB = 5 - 2/ - 2/e, 



a) 



/C(nB), /C(nB) 



97r 



(2.7) 



where (ki^/aH), as aheady stressed after Eq. (II. 3p . in the ACDM scenario, is solely deter- 
mined in terms of e and An- In Eq. (12. 7p /C(nB) varies very littltH in the range of physical 
rates 2 < / < 2.2; for values / > 2.2 (excluded from the present considerations) the energy 
density fluctuations induced by the magnetic fields will get larger than the adiabatic mode. 
The amplitude AB{ki^,T) depends on the thermal history through the factor (oi/a) which 
gets different values depending on the epoch at which the power spectrum is evaluated, for 
instance, at equality, 

/2fiR^y/' [H^/nMo' 
'UeAn) \ Mp\n 



ai 
a, 



eq 



RO 



(2i 



MO 



where /iq^ro = 4.15 x 10~^ and the total matter fraction at the present time is Q 
Qco + ^bo is given by the sum of the CDM and of the baryonic contribution. Following the 
same logic of Eq. (12. 8p . but at a different redshift, Eq. (II. 4p can be derived from Eq. (12. 7p 
by recalling known conversion factors between different system of unitJ^. 

The parameters of Eq. (12.80 can all be extracted, assuming the ACDM model, by using 
different data sets [UllTTlIIH] (see also e.g. Eqs. (14.30 - 04.50 in section Hj). The only exception 
is represented by e which cannot be determined in the vanilla ACDM model but it can be 
bounded from above in terms the tensor to scalar ratio rx = Vt/V-jz measuring the ratio 
between the spectrum of the tensor modes and the spectrum of the curvature perturbations 
at the pivot scale kp. To first-order in the slow- roll expansion the tensor to scalar ratio 
Tt = 16e + 0{e^) so that any bound on rx can be translated into a bound on e. The 



Data 




ris 


e 


WMAP9 alone 


< 0.38 


0.992 ±0.019 


< 0.023 


WMAP9 + Hubble 


< 0.34 


0.995 ±0.015 


< 0.021 


WMAP9+ BAG 


< 0.18 


0.973 ±0.011 


< 0.011 


WMAP9+ all 


< 0.13 


U.yU^( _Q QQg4 


< 0.0081 



Table 1: The upper limits on the tensor-to-scalar ratio and on the slow- roll parameter for 
some illustrative choices of cosmological data sets. 



WMAP9 observations can be combined with different data sets in the light of the ACDM 
scenario supplemented by the tensor modes of the geometry. In this way it is possible to 

^Note that /C(l) = 1. 

^"It is useful, for numerical estimates, to appreciate that Mp = 2.149 x 10^^ nG. 
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obtain a bound on and on e. In Tab. [T] the bound on rx obtainable in the case of the 
WMAP9 data alone is compared with the same bound obtained by combining the WMAP9 
data with the measurements on the Hubble constant of Ref. [13] or by combining the 
WMAP9 observations with the baryon acoustic oscillations, dubbed BAO in Tab. [1] (see e.g. 
[22]). The last row of Tab. [T] refers to the combination of the WMAP9 data with almost all 
the data sets recently available, namely, the data on the Hubble rate, the ones on the baryon 
acoustic oscillations supplemented by the ones of the Atacama Cosmology Telescope [44] , 
by the data of the south pole telescope [15] and by the three year sample of the supernova 
legacy survey [46] . 

For the standard thermal history with sudden reheating also the conductivity cXc jumps 
at a finite value at the end of inflation and the continuity of the electric and magnetic fields 
implies that the amplitude of the electric power spectrum gets suppressed, at a fixed time, 
as {k/acY comparison with its magnetic counterpart [35]. The electric mode functions 
are exponentially suppressed, for a fixed wavenumber. The magnetic power spectrum is also 
suppressed, for sufficiently large A;, as exp [— 2(A;^/A;^)] where /c^^ = (ir'/[47r(Tc(r')]. The 
evaluation of is complicated by the fact that the integral extends well after Tq-. This 
estimate can be made rather accurate by computing the transport coefficients of the plasma 
in different regimes [47]. By taking r = Tcq 

4.75x10-^^ / k X- 

showing that exp [—2{k/kaY] ^o small to give neghgible suppression for 0{kp) < k < 0{ki^) 
where the present considerations apply. 



2.2 Energy density and anisotropic stress 

The normalized fluctuation of the magnetic energy density and of the related anisotropic 
stress are given by: 

Spb{x,t) = J j0j-,SpB{q,T)e-^^-, nlf (f,r) = / ^^njf (g»e-^'?-, (2.10) 

The normalized fluctuation of the magnetic pressure is 6pb{x,t) = 6pb{x,t) /3. Since we 
shall be dealing with the scalar modes of the geometry, it is practical to introduce the scalar 
projections of the anisotropic stress 

V^UB{x,T) = d,d,Ii:^{x,T). (2.11) 

The fluctuations of the energy density and of the anisotropic stress in Fourier space are 
reported in appendix|X]to avoid lengthy digressions since this analysis is standard and can be 
found within slightly different perspectives in [191 [231 136]. In full analogy with the magnetic 
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fluctuations the fluctuations of the electric energy density and of the electric anisotropic 
stress can be defined. 

The magnetic energy density and and the anisotropic stress induce scalar fluctuations of 
the geometry as established long ago (see [19] and references therein). The variables h{k,T) 
and ^{k,T) parametrize the metric fluctuations in the synchronous coordinate system 

6, 



5sgij{k,T) = a 



kikjh -\- 6^1 ki kj 



Ss9oi{k, t) = Ssgoo{k, t) = 0, 



(2.12) 



where 6s stresses that we are here considering the scalar modes of the geometry. For practical 
reasons, in what follows, this subscript will be omitted. 

The evolution of h{k, r) and ^{k, r) can be obtained by perturbing the Einstein equations 
during the inflationary phase: 



-2k^i + m' = -ATTGa\6p^ + 6p^ + 6pb), 
h" + 2Hh' - 2k'^i = 2A7rGa\6p^ + 6p^ + dps) 
{h + 60" + 2n{h + 60' - 2k'^^ = 247rGa^nB, 
e^' = -ATxG[ip'k\^ + a'k\, - P], 



(2.13) 
(2.14) 
(2.15) 
(2.16) 



where x^p and Xa denote, respectively, the fluctuations of ip and of a. Equations f l2.13p 
and fl2.16p come, respectively, from the (00) and (Oz) components of the perturbed Einstein 
equations. Equations (I2.14p and fl2.15p are derived from {i = j) and [i 7^ j) components of 
the perturbed Einstein equations. In Eqs. (12.151) and (12.161) the three-divergence of both 
sides of the equations has been taken. Finally, in Eq. fl2.16D P{k, r) denotes the Fourier 
transform of the three-divergence of the Poynting vectoi^l], i.e. V ■ [-E x 5]/(4'7ra^). The 
electric fields have been neglected since they are suppressed all along the inflationary phase. 

In the synchronous gauge description, the fluctuations of the energy density and pressure 
of (f and a are: 



Sp^ 



1 




^ dV 






dip 


1 




dV 






dip 



d Xifi 



Spa 

Spa 



dx'. + 



/ / 

(^Xa 



dW 

dW 

dip 



(2.17) 
(2.18) 



Te fluctuations of the inflaton and of the spectator field obey, respectively, the following two 
equations: 



X'' + 2Hx!^ + ^'X^ + ^a'x^ -^h' = 



x: + 2nx'^ + k'x. 



dcp'^ 
d^W 
~d^ 



a^XcT 



-h + — 5pB 



0. 



(2.19) 
(2.20) 



^^While P can be neglected in comparison with the other contributions of the momentum constraint, 
P' = —AHP + k^{SpB — 3nB)/3. Even if P has been neglected in the final expressions, its derivative has 
been used whenever needed. 
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2.3 Curvature perturbations 

The system of the perturbed Einstein equations f l2.13p - fl2.16p supplemented by Eqs. f l2.19p 
and fl2.20p can be decoupled in terms of two variables defined as: 

= aXip - z^i^ , g<7 = axa - z^i, (2.21) 

where z^p = aip'/H and z^- = acr'/'H. The evolution equation obeyed by and are 
obtained from Eqs. (KWl and (^M> with the help of Eqs. and fl^TTB]) . The general 

result can be further simplified under the assumption that the energy density of a can be 
neglected in comparison with the energy density of the inflaton. Defining the rescaled Planck 
mas^ the Friedmann equations 

SMln^ = ^{ip'^ + a'^) + a^V{if) + a^W{a), 2mI {H^ - H') = ^'^ + a'\ (2.22) 
imply that the evolution equations for q^ and g^- can be expressed as 

% + k% -^q^ + ^ S^{k, r) = 0, (2.23) 

q'; + k'q. - + ^ :s.(fc, r) = 0. (2.24) 
a M-p 

Equations f l2.23p and f l2.24p hold under the approximation that z^ <^ z^ (as implied by the 
occurrence that a' <^ (f'); the source terms in Eqs. f l2.23p and f l2.24p are given by 

The curvature perturbations on comoving orthogonal hypersurfaces are expressible, in 
the synchronous gauge, solely in terms of ^ and of its first time derivative, i.e. 



Using then Eq. fl2:T6|) into Eq. (jMID to ehminate f and recalling Eq. ([MI]), 7^(A;, r) 
becomes: 

V(kT]- ^^i'^)%ik,r) + z„{T)q,{k,r) _ q^{k,r) .z,{r) , . 

where the second equality follows again in the limit z^ z^p. Equations fl2.25p and f l2.26p are 
easily solvable with standard Green's functions methods when the relevant wavelengths are 



-"^^In the present paper we shall use both Mp and Afp. The two quantities are equal up to the factor v^Stt, 
i.e. Mp = Mp/Vs^ I/VS^. 
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larger than the Hubble radius during inflation; the result of this procedure can be written 



as 



qaik, t) 



c^flB{k, t) + d^VL-Bn{k, r) 
c^^Bik, r) + d^QBu{k, r) 



a(r), 



(2.29) 
(2.30) 



where the sources have been evaluated to leading order in kr and 



^B{k,T) 



6pB{k,T) 



f2 ' 



^Bn{k,T) 



llB{k,T) 



In Eqs. fl2.29p and ( I2.30p the following quantities have been introduced 

3 



(2.31) 



3Mp \HJ' 
1 



m{f,e) 



d = -c 



dn 



2Mp 



m(/,e) 



3(1 



(l-2/)(4-2/-3e)- 



(2.32) 



(2.33) 



Using Eq. f l2.29p and f l2.30p inside Eq. f l2.28p the resulting expression for the curvature 
perturbations can be written as: 



n{k,T) = n,{k) + S^ik) + MfXr)^B{k,T) + M^^J'\T)QBu{k,T 

where TZ^{k) denotes the standard adiabatic solution associated with q^\k^T 
non-adiabatic mode associated with q^\k,T): 



n,{k) 



S.{k) = -ql'\k,T) 



(2.34) 
S^{k) is the 

(2.35) 



The two functions A^^^^(r) and M.^?'^\t) are defined as 



Mpa{T) 



Mp air) 



d^ + 



dn 



(2.36) 



The result obtained in Eq. fl2.34p has been deduced in rather general terms and it holds, with 
different forms of the coefficients, not only when the gauge fields are coupled to the spectator 
field but also when some direct coupling to the inflaton is present. In the forthcoming sections 
we shall assume the result of Eq. fl2.34p and keep the coefficients general. In terms of these 
coefficients the initial conditions of the Einstein-Boltzmann hierarchy can be derived. In the 
specific case discussed here the dominant contributions will be the ones of c^p and d^ which 
are, at most, 0{l/e) in the nearly scale-invariant limit. 
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2.4 Global variables 



There is a conservation law associated with the evolution of TZ. Recalling the explicit ex- 
pression of Eq. fl2.27p . we can take the first derivative of both sides and use the evolution 
equations of the fluctuations fl2.13l) - fl2.16l) where, instead of specifying the energy densi- 
ties and the pressures we simply introduce 6pt (the fluctuation of the total energy density), 
6pt (the fluctuation of the total pressure) and the total non-adiabatic pressure fluctuation 
5pnad = Spt — c^t^pt where c^^ = p[/ p[ is the total sound speed of the system; the background 
pressure and energy density pt and pt obey the conventional Friedmann equations in the 
spatially flat case: 

mlH^ = a^Pt, 2mI{H^ -W) = a\pt+pt). (2.37) 
Equation fl2.37p implies p[ + 3'H(pt + Pt) = 0. The first derivative of TZ turns out to be 

^ _mpnad ^ n5pB / 2^ _ 1 \ n\h' + 6^) _ uk^cU ^ uii, 

Pt+Pt Pt + PtV''* 3/ 87rGa2(pt + Pt) AnGa'^ipt + pt) {Pt + pt)' 

where lit denote the total anisotropic stress of the system containing together with the 
anisotropic stress of the magnetic fields also the anisotropic stress of the other species. In 
the ACDM scenario the other source of anisotropic stress comes from the neutrino sector. 

If S^{k) = in Eq. f l2.34p . 5pnad = in Eq. (12. 38 p . This will be the situation considered 
in the forthcoming section when discussing the magnetized adiabatic mode. Recalling the 
results of the first paper of Ref. [23], the presence of the non-adiabatic modes can also 
be easily considered along the lines illustrated above but for the purposes of the present 
discussion it is unnecessary. 

By taking the time derivative of both sides of Eq. (I2.38P and from the evolution equations 
for h and C, the equation of TZ simplifies even further: 

7^" + 2^7^' - 4 V'7^ = + 2^S7^ + ^Ht, (2.39) 

Zt Zt 



where 



— ~- : r + 



{Pt + Pt) (pt + Pt) 



4 - ^ )^pb + nt 



= -y^, (2.40) 



The results of Eqs. (I2.39p - (l2.40p reduce, in the absence of magnetic fields and in the absence 
of anisotropic stress, to the results of Lukash valid in the case of an irrotational relativistic 
fluid pi| |2ni For wavelengths larger than the Hubble radius the Laplacian of TZ can be 
neglected and the solution of Eq. fl2.39p is 

7^(f,r) =7^,(f)+ / dT'T.n{x,T')+ a"(r') nt(f , r') rfr. (2.41) 
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The continuity of TZ across the inflation-radiation transition can be verified in exphcit toy 
models. Consider, for instance, the situation where pt and pt are both continuous across the 
inflation-radiation transition. Denoting with /3 = a/a* the normalized scale factor across 
the transition, we must require that the effective barotropic index Wt — )■ —1 and Wt — )• 1/3 
when /3 gets, respectively, much smaller and much larger than 1. An interesting interpolating 
solution of Friedmann equations with this property is given by 



_ 12H'^ Ml _ 2 iP' - 3) 

A - (^^^' - 4^*^P (^^- (2-42) 

The barotropic index Wt = {l/3){l3^ -3)/{l3'^ + 1) goes to 1/3 for /3 > 1 and to -1 for /3 < 1. 
In conformal time the evolution equations can be explicitly integrated with the result that 
H = (r^ + T^)"^^^ and /3(r) = + ^ + . These expressions can be used to verify the 

continuity of 7l{x, r) for instance, in Eq. (12.41 p . 

When inflationary magnetic fields and non-adiabatic pressure fluctuations are both van- 
ishing, the initial conditions of the temperature and polarization anisotropies follow from 
Eqs. f l2.39p . f l2.40p and fl2.4ip . Even if inflationary magnetic fields are absent the total 
anisotropic stress lit receives contribution from the neutrinos. The compatibility of the neu- 
trino anisotropic stress with the other evolution equations implies, in the case of adiabatic 
initial conditions, the well known result [10| HI] stipulating that lit — 0{k'^T^) when the rel- 
evant wavelengths are larger than the Hubble radius at the corresponding epoch. Using Eq. 
f l2.4ip and the interpolating solution of Eq. f l2.42p it can be shown by direct integration that 
during the radiation epoch 7^(A;, r) ~ 7^*(A;){1 - (4/9) [/?,./ (4/2^ + 15)]Pt'^} where Ru < I 
is the (massless) neutrino fraction of the radiation plasma. In summary, for adiabatic initial 
conditions and it the absence of magnetic fields of inflationary origin Tl{k, r) ~ TZ^{k) with 
corrections which are always small for typical wavelengths larger than the Hubble radius. 
The generalization of this statement to the case when the inflationary seeds are present will 
be, among other things, the subject of the following section. 



3 Initial conditions of the CMB anisotropies 

The value of curvature perturbations computed in the previous section will now be used as 
initial condition for the subsequent post-inflationary evolution. Focussing on the magnetized 
adiabatic mode the post-inflationary value of curvature perturbations will be parametrized 
as: 

7^(A;, r) = n^{k) + ceR^^Bik) + dBRy(rB{k), (3.1) 

where cb and ds are of the same order and, as discussed before, at most 0{l/e). In Eq. (13. ip 
the energy density and the anisotropic stress of the magnetic fields have been rescaled through 
the energy density of the photons. After neutrino decoupling (occurring approximately 
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around the MeV) the anisotropic stress of the neutrinos is generated and Ilt{k,T) can be 
parametrized, in the ACDM scenario, as 

Ilt{k, t) = {p^ + py)au{k, t) + {p^ + p^)crB(A;), (3.2) 

where, following the conventions of section [21 a^, is related to the anisotropic stress in real 
space as didjWJ = {p^ + p^) V^cXj,. 

3.1 Pre-decoupling plasma 

The perturbed Einstein equations for the post-inflationary system of magnetized perturba- 
tions include photons, neutrinos, baryona^i and CDM particles. The analog of Eqs. f l2.13p - 
fl216|) will then be: 

2e^-nh' = £^{6pt + 6pB), (3.3) 
Mp 

3fl2 

h" + 2'Hh' -2e^= ^{Spt + SpB), (3.4) 
Mp 

{h + 60" + 2n{h + 60' - 2k^^ = =^[ipu + Pu)cTu + (p^ + p^)aB], (3.5) 



a 



2 



k'C' = -^{Pt + Pt)0t, (3.6) 
2Mp 

where 9t, 6pt and 6pt ar 

(Pt + Pt)Ot = + Pa)6'a, Spt = 5sPa, ^Pt = '^sPa = ^a^sPa- (3.7) 

The sums appearing in Eq. (13. 7p extend over the four species of the plasma (i.e. photons, 
neutrinos, baryons and CDM particles) and Wa is the barotropic index of each species. By 
using the background equations in their general form, i.e. Eq. f l2.37p . Eqs. (13. 3p and (13. 6p 
can also be written in more explicit terms as: 

2k^^ - m' = Sn^hniR-ySy + RJ.) + R^^r^b + ^^iviif^Vc + (77^ Vbl |(3-8) 

L LViZMO/ V"MO/ J J 



3 / ^ 

^kRv(^v + ^R-R7(l + -Rb)^7b + T^M 7:; ]PcOc 

4 ViiMO- 



(3.9) 



where 5^, 5^, (5b and 5^ denote, with obvious notations, the density contrasts of the four 
species of the plasma. 



-"^^The evolution of electrons and ions can be described, effectively, in terms of a single fluid called the 
baryon fluid whose velocity is the center of mass velocity of the electron-ion system. This happens because 
the Coulomb coupling is always tight, as argued long ago in the pioneering work of Peebles and Yu |39| . 

^'''In Fourier space 9t and 9^ denote, respectively, the three-divergences of the total velocity field and of 
the different species composing the plasma. 
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At sufficiently early times, the velocity of the photons coincides with the baryon velocity, 
i.e. 6^^ ~ 6'b = since the photons and the baryons are tightly coupled by Thomson 
scattering. In the momentum constraint of Eq. (13. 9 p the tight- coupling regime has been 
already assumed. Furthermore, because we are going to solve the system across the radiation- 
matter transition, the Poynting vector can be safely neglected. Using the same strategy 
leading to Eqs. (13. 8 p and (13. 9p . Eqs. (13. 4p and (13. 5 p can be recast in the following form: 



{h + 60" + m^h + 60' - 2fc2^ = l2H^n^[Ru(yy + R^^b] 



(3.10) 
(3.11) 



where, the critical fraction of radiation and of matter are defined as: 



n 



R 



a 



M 



a + a, 



(3.12) 



eq 



The initial conditions of the Einstein-Boltzmann hierarchy are fully specified by solving 
also the equations of the monopoles and the dipoles of the phase-space distribution of the 
various species. The monopoles are related to the density contrasts whose evolution is given 
by: 



5', 



4 2 
-9u + -h', 
3 3 ' 



4 2 

3 ^ 3 ' 



77b 



h' 
~2' 



SI = + 



h' 



(3.13) 
(3.14) 



The equations for the dipoles are related to the peculiar velocity of the various species: 



e:. 



9' + Wc = 0. 



-5^ + 



4(l + i?b) ^ 4(l + i?b) 



(^B - 4aB) 



(3.15) 

(3.16) 
(3.17) 



In Eq. (I3.15P there appears also a^, whose evolution equation is coupled to all the higher 
multipoles of the neutrino phase-space distribution 



4 3 2 4 

- —kT,^ - —h' - -0, 

15 10 15 5^ ' 



k 



vl 



21+1 



(J^iy{e-i) — (^ + ^)J^u{e+i 



> 3. 



(3.18) 
(3.19) 



where J-'^i denotes the £th multipole of the perturbed phase-space distribution of the neu- 
trinos. The evolution equations reported here hold prior to photon decoupling. Across 
decoupling the baryons and the photons obey effectively different equations since the ap- 
proximation based on the tight photon-baryon coupling breaks down. 
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3.2 Explicit solutions and magnetized adiabatic mode 

The solution of Eq. fl2.37p across the radiation-matter transition and for a spatially flat 
Universe reads: 



« = - = x^ + 2a;, x = ^, -1 = ^\/|^- 283.73 (^)"mpc, (3.20) 

Oeq Ti iio V "MO V0.1368/ 

where Oeq is the scale factor at the equality already introduced in Eq. fl3.12p . i.e. the moment 
when non-relativistic matter and radiation contribute equally to the total energy density of 
the plasma. For a = pu/ Pr ^ 1 (i-e. a ^ a^q) the plasma is dominated by radiation and 
according to Eq. ^Ml, a 2x + 0{x^) = 2{t/ti). 

Defining as Ti the initial integration time, it will be required that kr-i < 1 for all the modes 
involved in the calculations. The double expansion employed in setting initial conditions of 
the truncated Einstein-Boltzmann hierarchy can be formally expressed aJ^ 

k k Ki (y. ^ , , 

aH ti 2 V a + 1 

where k = kri measures how large the wavelength was, in Hubble units, around equality 
(note, indeed, that r^q = {\/2 — l)ri ti/2). 

The initial conditions studied here incorporate the inflationary seeds in the adiabatic 
mode and belong to the class of the magnetized adiabatic modes (see [HI [23] and references 
therein). To investigate the time evolution of the system it is useful to employ directly the 
normalized scale factor a. The initial conditions pertaining to the magnetized adiabatic 
mode can then be written as: 

3 3 

S^{k, fti) ~ S^{k, fti) ~ -Shin, Oi) ^ -5c = ~R-y^Bifi, tti), (3.22) 
9u{k, ai) ~ Ojhifi, tti) - Oc{fi, tti) - 0, (3.23) 
a^{K,ai) = 0, J^i{K,,ai) = 0, (3.24) 

where ai denotes the initial integration variable and i > 3. In the synchronous gauge, the 
density contrasts on uniform curvature hypersurfaces are defined as Ca = + '^a/[3(wa + 1)]- 
But then, the conditions f l3.22p imply that all all the Ca niust be equal, i.e. Ci/ = C7 = Cc = Cb 
for a = a-i and k < 1. In this sense the initial conditions fl3.22p- f l3.24p generalize the 
adiabatic mode to the situation where inflationary seeds are pre33. 

According to Eq. fl2.39p the initial conditions f l3.22p - fl3.24p seem to be compatible with 
the adiabatic mode only if lit — > 0. The neutrino anisotropic stress must be zero initially 
and also its flrst derivative is zero since J-'^e = for i > 3. But thanks to the presence 



^^In terms of the solution of Eq. (|3?20| . H ^ 2y/a + l/{Tia) and c^^ = A/{3a + 4). 

^^It should be borne in mind that, in the synchronous gauge, the condition 0c = is enforced not so much 
because of a property of the initial data but rather to fix completely the coordinate system and to avoid the 
occurrence of known spurious (gauge) modes arising in the synchronous description [42] 
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of the magnetic anisotropic stress, the total anisotropic stress does not vanish even before 
neutrino decouphng: Ilt{K,ai) 7^ even if (T,y(K, Oi) and its derivatives are all vanishing. 
Direct numerical integration shows that, after a transient time, lit — even if, initially, 
Ut ^ 0. This result is established by integrating Eqs. (13:41) - fim]) . (I37[3|) - (l3:ii|) and f l3J[5|) - 
f l3.19p in the background defined by Eq. fl3.20p and subjected to the initial conditions (13.1 p 
and fi3.22p - fi3.24p . In Fig. [T] the evolution of the neutrino anisotropic (T^{K,a) stress and 



a „=0.0463, £2 „=0.233, £1^ =0.721 
bO cO A 



=0.0463, a „ =0.233, =0.721 

bO cO A 




K =0.001, 0.01, 0.1 



7 -6 -5 -4 -3 -2 -1 




loga 



log a 



Figure 1: Evolution of the anisotropic stress in the case of the magnetized adiabatic mode. 
In the left plot the neutrino anisotropic stress is illustrated. In the plot at the right the total 
anisotropic stress is reported. On the horizontal axes the results are given in terms of the 
common logarithm of the normalized scale factor. 

of the total anisotropic stress Ilt{K,,a) is illustrated as a function of the common logarithm 
of the normalized scale factor. On the vertical axis of both plots of Fig. [1] the anisotropic 
stresses are given in units of the magnetic anisotropic stress o"b(k). 



Denoting with R,, and 



Ri, the fractions of neutrinos and photons in the pre- 



decoupling plasma^l. the arrow in the left plot underlines the value —R^/Ry (which is asymp- 
tote of ay/a-Q for kr <^1). The result of Fig. [T]can be obtained by integrating directly the 
system in the a parametrization. The equations are reported, for completeness, in appendix 
[Bl As the legends indicate in Fig. (H different values of n produce results which are indistin- 
guishable provided k <^\. 

In the limit a < 1 and A;r < 1 we can solve consistently the system discussed in the 
previous section and expressed, in the a-parametrization, in appendix [Bl The result is: 



12(4i?, + 15) 



1 - 



+ R^n^{k) 



Ryk^r' 



cb 



6{ARy + 15) 



i^In the ACDM paradigm i?^ = [3 x (7/8) x (4/ll)4/3]/[l + 3 x (7/8) x (4/11)4/3] ^ 0.4052 where 3 counts 
the massless neutrino famihes, (7/8) stems from the Fermi-Dirac statistics and (4/11)''/^ jg related to the 
different kinetic temperature of neutrinos. 
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h{k,r) = 

+ 



dB + 



TZ:t,[k) 2 2 

K T 



2Pr 



2^2 



3(4i?, + 15). 
^ ^ 8i?2 - UR^ - 75 



36(2i?^ + 25)(4i?^ + 15) 
rJ{20R^ - 15) 



Cb + 



180(4i?^ + 15)(2i?^ + 25) 



^ i^O-B(fe) ^^2^2 



i?^(20i?^ - 15) 



45(4i?, + 15)(2i?, + 25) 



For the density contrasts we have instead 



T^*{k) , 2 2 



k'T'-R^n^{k) 



1-77 



2R^ 



k^T^ 



1 + 



dsR^ 



2_2 



5^(A;,t) = 



TZ^ik) 2 9 



3 



l + {R^-2c 



B 



df 



2R, 



2_2 



k'r 



Sc{k,T) 



+ 



4 4 ' ^ ^ 

^k'r\2 + dsR,) 



k^T^ 



1 - Cb- 



+ 



dsR-y 



aB{k)k'r' 



2i?, 



For the velocities we have 



9-yh{k,T) 



36 



9V i?. 



1 + ^(CB^-^)^V2 



- CTBik) 



1 + 



1 



1 - 



df^R^ 



2 

4^3 



\ ^2^2 



(3.25) 



(3.26) 



(3.27) 

(3.28) 
(3.29) 

(3.30) 



i^v > y *v ^V4i?^ + 15y 36 4 



^ _ 1 /4i?, + 23 



9Ui?^ + 15^-^ 2(4i?^ + 15) 



4i?, + 27 



+ C7B(k)^k\\l 



k^T^ 



AR^ + 27 deR^ f4:R^ + 23 



2(4i?, + 15) 



+ 



AR^ + 15 



(3.31) 
(3.32) 



while 9c = 0. Finally the anisotropic stress of the neutrinos is given by 
(Tj.(A;,r) = 



r/^^"^ 3(4i?. + 15)^^ 



3(4i?, + 15) 



CB + -]R^Vt^{k) 



2a^{k) R^ / 



A';V'(3.33) 



(4i?, + 15) i?. 

All the higher multipoles in the neutrino hierarchy have been consistently set to zero since we 
are concerned here with the adiabatic initial data. Different initial conditions may demand 
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different assumptions on the higher multipoles of the hierarchy. This analysis closely follows 
earlier results (see [191 123] and references therein) where, for the first time, the problem of 
the magnetized initial conditions of the Einstein-Boltzmann hierarchy has been posed and 
discussed. The crucial difference is represented by the inflationary origin of the magnetic 
fields which induces a further contribution on the curvature perturbations and, therefore, on 
the whole hierarchy. Note that when the magnetized contribution is switched off, the initial 
conditions of Eqs. fl3.25p - fl3.33p reproduce the standard adiabatic initial condition discussed 
long ago (see, e.g. |^). 



4 Temperature and polarization power spectra 

The temperature and polarization observables can be obtained by integrating numerically 
the magnetized Einstein-Boltzmann hierarchy across recombination with the full set of initial 
conditions discussed from Eq. f l3.25p to Eq. 03.33 p . The Boltzmann integrator employed 
here is based on the code described in Ref. [19], [23] and used to investigate the magnetized 
CMB anisotropics. The Boltzmann code is based on Cosmics [lOj and on CMBFAST [H] 
and it includes the evolution of magnetic fields within a consistent magnetohydrodynamical 
approximation. 

The temperature autocorrelations (TT correlations in what follows), the polarization 
autocorrelations (EE correlations in what follows) and the cross- correlation between the 
temperature and the polarization (TE correlations in what follows) are defined in the stan- 
dard standard way (see, for instance, Eqs. (2.61)-(2.63) of the second paper quoted in Ref. 
|23]). All the different correlation spectra will be expressed in units of (/uK)^. The following 
shorthand notation shall be used: 

(TT) _ £(£+1) TT ^(EE) _ ^(^+1) ^EE ^(TE) _ ^(^+1) ^TE ^ 1 ^ 

/TT /Vr ZTT 

In the minimal ACDM scenario the angular power spectra of Eq. (14. ip are the only non- 
vanishing observables since the tensor modes are neglected and the B-mode polarization is 
absent. As we shall argue later, however, a B-mode autocorrelation can be induced, via 
Faraday effect, from the EE correlations. 

The spectra of Eq. (14. ip depend on 6 independent parameters 

Qf'^^ = Qf^^\ns, Qho, ^co,^A, Ho, ere), (4.2) 

where ere (not to be confused with the slow-roll parameter) denotes the optical depth at 
reionization. The pivotal parameters of the ACDM paradigm can be determined on the 
basis of different data sets and, for illustrative purposes, we shall bound the attention only 
to three best fits. The first one is obtained by comparing the ACDM paradigm to the WMAP 
9yr data alone (see, in particular, [1^): 

(fibo, f^cc^^deo, ho, ris, ere) = (0.0463, 0.233, 0.721, 0.700, 0.972, 0.089), (4.3) 
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a„ =0.0463, n „ =0.233, = 0.721 a „ = 0.0463, a „ =0.233, =0.721 

bO cO A bO cO A 




I I 

Figure 2: In the two upper plots the TT and the EE correlations are illustrated for the initial 
conditions of Eqs. f l33]) -( l333|) 

with A-ji = 2.41 X 10~^ (recall, in fact, the parametrization introduced in Eq. fll.2p ). If 
we include the data sets pertaining to the baryon acoustic oscillations (see, e.g. ^2]) the 
parameters are slightly different: 

(fibo, fico,^deo, ho, n^, e^e) = (0.0477, 0.247, 0.705, 0.686, 0.967, 0.086), (4.4) 

with Atz = 2.35 x 10^^. Another possible set of parameters considered hereunder is the one 
obtained by combining the WMAP9 data with the direct determinations of the Hubble rate 

(fibo, f^cc^^deo, ho, n„ ere) = (0.0445, 0.216, 0.740, 0.717, 0.980, 0.092), (4.5) 

with An = 2.45 x lO'^. 

The presence of inflationary magnetic fields affects the CMB observables obtained in the 
framework of a particular best fit to the WMAP 9yr data for a sufficiently large value of e. 
This aspect is illustrated in Figs. |2] and |3] where the magnetic fields have been introduced both 
at the level of the initial conditions and at the level of the evolution equations, as discussed 
in section |3l Both in Figs. |2] and |3] the fiducial set of ACDM parameters has been chosen 
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a „ =0.0463, n = 0.233, =0.721 a „ =0.0463, Q. = 0.233, =0.721 
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Figure 3: The same cases of Fig. [2] are illustrated in terms of the cross-correlation between 
temperature and E-mode polarization. At the right the detail of the first anticorrelation 
peak is illustrated. 

as in Eq. (14. 3p . In the framework of the ACDM scenario with standard thermal history 
the only extra parameter is represented by the slow-roll parameter e: the other parameters 
controlling the amplitude of the magnetic power spectrum of Eq. (12. 7p are determined by 
the other parameters of the concordance scenario. 

In Fig. [2] for sake of simplicity two extreme examples are illustrated, i.e. the case 
e = 0.04 (full line) and e = 0.001 (dashed line). In the two upper plots of Fig. |2] the 
TT and EE power spectra are reported, while in the two lower plots of Fig. [2] the angular 
power spectra appearing in the upper plots have been divided by the best fit to the vanilla 
ACDM model (i.e. in the absence of inflationary magnetic fields). Using the notations of Eq. 
(14. ip in the two lower plots of Fig. [2] we illustrated, respectively, Qi^^\e) / Q^^^\e = 0) and 
^£™''(e)/^£™''(e = 0). The shorthand notation e = simply means that the corresponding 
power spectrum is taken to be independent of e, as it happens in the vanilla ACDM with no 
tensors. In Fig. |3] the same analysis has been performed in the case of the TE correlations 
which cannot be simply divided by the corresponding WMAP 9yr best fit since the TE 
correlations are not positive definite. To illustrate more closely the differences between the 
different models, in the right plot of Fig. |3] the first anticorrelation peak has been shown in 
greater detail. 

From Figs. |2] and [3] the case e = 0.001 is practically indistinguishable from the WMAP 
9yr bestfit while the case e = 0.04 shows quantitive and qualitative differences potentially 
jeopardizing the agreement of the computed spectra with the observational data. This 
disagreement arises since the first, second and third peaks of the acoustic oscillations are 
distorted. This aspect can be scrutinized from the values of the TT correlations in the 
neighborhood of the first three acoustic peaks. The position of the peaks will be denoted, 
respectively, by ii, £2 and £3. To pin down the position of the peaks in terms of the ACDM 
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parameters we use the following parametrization adapted to WMAP 9yr data |23] : 



£,=e,+Ai,, = £^(j-_0.), (4.6) 

where (pj and Aij are given, for j = 1, 2, 3, as follows: 

0, = 0.267 (^)° \ 02 = 0.241 (^)° ', 03 = 0.353 (^)°'', (4.7) 
Ail = 0.13 \n,-l\Ii, A£2 = 0.33|n, -1|£2, Ai^ = 0.61\n, - 1%. (4.8) 

Note that is simply the well known acoustic multipole which is expressible through the 
angular diameter distance to recombination. Its standard expression can be reduced to a 
more explicit formula 



In 



v^T+J?^+v/(l+rR,)i?.b, 



i+Vm*i?b. 



(4.9) 



where tr.,, and -Rb* are given by 



TR. = = 4-15 X \0-\hl^^)-' f-^), R^{z) = ^^=30.36hln,o (—). (4.10) 

Pm[z*) VIO'^/ 4p^ V 2* / 

The quantity z^, is the redshift to recombination which can be directly expressed in terms of 
ACDM parameters as 

= 1048[l + (1.24x 10-=^)(/i2fibo)-°-'''][l+^?i(/i^fiMo) ^1, (4.11) 
0.0783 (/igr^b)""-^'' 0.560 
[1 + 39.5 (hln^o) ' 1 + 21.1 (/i^n^^^) 1.81 • ^ ■ ) 

The parameters of Eq. (14.31) imply = 1090.95 in excellent agreement with the estimate of 
Ref. [13 [IS] (i.e. z^ = 1090.41 ± 0.57) in the case of the WMAP 9yr data alone in the 
light of the vanilla ACDM scenario. The relative heights of the acoustic peaks computed in 
the case of the best-fit of Eq. (14. 3 p are: 



ifi = = 6.942, 1/2 = ^^ = 0.447, H3 = ^^ = 0.981, (4.13) 



where £i = 221, £2 = 538 and £3 = 815 are, respectively, the locations of the first three 
acoustic peaks obtainable from Eqs. (14. 6 p - (14. 91) and coincide, approximately, with the ones 
directly obtainable from the angular power spectra. In the case of the angular power spectra 
illustrated with the full line in the left plots of Fig. [2] the same ratios of Eq. (I4.13P are 
Hi = 7.35, H2 = 0.453 and H3 = 0.976. This shows that the largest effect is on the first 
peak while the other two are comparatively less affected. 

The results obtained in the case of the parameters of Eq. (14. 3 p are quantitatively and 
qualitatively close to the results obtainable in the cases of Eqs. ( 14. 4p or ( 14. 5p . In Fig. [H 
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Figure 4: The temperature and polarization autocorrelations divided by the corresponding 
best fits are illustrated for the fiducial set of parameters reported in Eq. fl4.5p . 

for instance, the parameters of Eq. (14. 5 p have been assumed while one of the values of e is 
different. From a closer comparison of Figs. [2] and HI it can be argued that for e < 0.03 the 
difference between the magnetized angular power and the corresponding ACDM best fit is 
smaller than 10^^. 

In the ACDM paradigm without tensors there are no sources of B-mode polarization but 
a stochastic magnetic field can rotate the polarization plane of the CMB at a rate depending 
on the difference between the refractive indices associated, respectively, with the positive and 
with the negative helicities. Faraday rotation is one of the situations where the inadequacy 
of the one- fluid approximation (for the baryon-lepton fluid) is manifest. The positive and 
negative helicities composing the (linear) CMB polarization experience, in a background 
magnetic field, two different phase velocities, two different dielectric contants and, ultimately, 
two different refractive indices. The mismatch between the refractive index of the positive 
and negative helicities induces, effectively, a rotation of the CMB polarization and, hence, 
a B-mode. The inclusion of the Faraday effect in the treatment implies, physically, that the 
proton-electron fluid (sometimes dubbed as baryon fluid) should be treated as effectively 
composed by two different species, i.e. the electrons and the ions. 

The Faraday effect for the CMB polarization can be treated either with uniform magnetic 
fields or with stochastic magnetic fields and different analyses have been performed starting 
with the ones of Ref. [H] (see also ^9] for an incomplete list of references). Since the 
inflationary magnetic fields does not break spatial isotropy, the angular power spectrum of 
Faraday rotation can be written as: 

{I{nMn,)) = i- j:{2i + l)Cf^P,{n, ■ h,\ X{n) = (4-14) 

where X{n) is a normalized form of the Faraday rotation rate, v denotes the comoving 
frequency and Pi{z) are the standard Legendre polynomials. In terms of the power spectrum 
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ft „ = 0.0463, a „ = 0.233, n =0.721 
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Figure 5: The B-mode polarization power spectrum induced by the Faraday rotation of the 
E-mode polarization. 

of the Faraday rate the autocorrelation of the B-mode polarization can be computed as 

E Z{i,iui2)Cf^Ci^^ (4.15) 

where Z{i,ii,i2) is a cornplicated function of the multipole moments containing also a 
Clebsch-Gordon coefficiento while Cf^ is the E-mode power spectrum already discussed 
above. 

In Fig. [5] the B-mode autocorrelation is computed in the case of the fiducial set of pa- 
rameters of Eq. (14.31) and for two illustrative values of e. The Faraday rotation rate depends 
also on the frequency of observation which has been taken to be 100 GHz. We remind that 
WMAP experiment observes the microwave sky in five frequency channels (i.e. 23, 33, 41, 61 
and 94 in units of GHz). The bandwidth increases from small to high frequencies signaling 
that probably the best sensitivity to polarization comes the high frequency channels. The 
WMAP 9yr data do not report any direct detection of the B-mode polarization. 

The Planck experiment |50j is observing the microwave sky in nine frequency channels: 
three frequency channels (i.e. u = 30, 44, 70 GHz) belong to the low frequency instrument 
(LFI); six channels (i.e. u = 100, 143, 217, 353, 545, 857 GHz) belong to the high frequency 
instrument (HFI). The BB power spectra for all the relevant frequency channels. There are 
reasons to expect that the sensitivity to polarization will be larger at high frequency [50]. 
At the same time the expected signal will be larger at small frequencies. We will illustrate 
the present results for a putative frequency of 100 GHz. 

-^^For an explicit expression of Z{£,£i,i2) see, for instance, the discussion contained in Appendix C of the 
last paper quoted in Ref. [23]. 
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The B-mode polarization induced by the tensor modes of the geometry is independent 
on the frequency channel and for rx ~ 0.1 the induced B-mode polarization ranges between 
10~^ (a^K)^ and 10~^(/iK)^ which is larger than what we have in Fig. |5]even assuming the 
maximum value of e. We must however bear in mind that the Faraday rate goes as z/~^ and, 
therefore, Cf^ goes as z/~^. This means that by moving from 100 GHz down to 30 GHz the 
signal roughly increases by a factor (10/3)'^ ~ 123. 

5 Concluding remarks 

Often, in the literature, sharp and challenging statements on the interplay of inflationary 
magnetic fields and CMB anisotropies are not corroborated by detailed analyses. We tried 
to explore here the opposite perspective by following a more pedantic approach involving 
the initial conditions of the Einstein-Boltzmann hierarchy in the presence of inflationary 
seeds. According to some, the really important problem is to tailor specific magnetogenesis 
mechanisms with very little attention to potentially interesting observational consequences. 
Some others prefer to assume the existence of large-scale magnetic fields and to scrutinize 
their phenomenological signatures without bothering to ask why those fields exist. As it has 
been shown in this paper, the two aforementioned perspectives must be seriously considered 
as complementary even if they are sometimes viewed as mutually exclusive. 

The necessity of bridging the gap between top-down and bottom-up approaches is, fortu- 
nately or unfortunately, inherent in all the analyses addressing the interplay of gravitational 
and gauge interactions in the early Universe. To be fair the same kind of problems emerge in 
the concordance lore. In this context concordance simply means a minimalistic agreement on 
the standard completion of the ACDM paradigm where the reheating is assumed to be almost 
sudden, inflation is driven by a single field and the post-inflationary history does not include 
long phases very different from radiation. Model independent analyses of the large-scale data 
are desirable but virtually impossible. In this spirit the initial conditions of the Einstein- 
Boltzmann hierarchy have been derived under the hypothesis that large-scale magnetic fields 
of inflationary origin were present prior to equality and for typical wavelengths larger than 
the Hubble radius at the corresponding epoch. Following earlier attempts, the guiding goal 
of the present investigation is to bring the treatment of magnetized CMB anisotropies to the 
same standards which are typical of the cases where large-scale magnetic fields are absent. 

The main assumption has been that the ACDM scenario with standard thermal history 
and inflationary completion is a sound approximation to the early dynamics of our Hubble 
patch. Here the amplitude of the magnetic fields is not given extrinsically as a further pa- 
rameter but it depends on the standard ACDM parameters. The parameters of the magnetic 
power spectrum and hence of the CMB observables only depend on the slow-roll parameter. 
As an application we presented the explicit computations of the temperature autocorre- 
lations, of the polarization autocorrelations and of the cross-correlation power spectra of 
temperature and polarization. B-mode autocorrelations are potentially generated by Fara- 
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day rotation of the CMB whose hnear polarization is affected, in turn, by the presence of 
the magnetic fields. 

Let us conclude with a conjecture. The ACDM scenario with tensor completion leads to 
an upper limit on the tensor to scalar ratio tt- Such a limit, by appropriately combining 
various data sets ranges from tt < 0.3 down to tt < 0.1. The limit on tt can be easily 
translated in a limit on the slow-roll parameter e and then we discover for instance that 
if tt < 0.12 the modifications induced by the inflationary seed on the temperature and 
polarization power spectra are indistinguishable from the differences associated with the use 
of different data bases. Conversely, if we believe that relic magnetic fields are there (and 
tensors are absent), then the present results suggest that e = 0.03 would already induce 
observable differences in the CMB spectra and this will imply an independent bound on e 
possibly achievable with more accurate analyses. In the near future less conventional models 
of inflationary magnetogenesis can be analyzed by using the same approach developed here. 
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A Correlation functions 



The comoving electric and magnetic fields are defined as 



E = c?^e, B = G?\f\h. 



(A.i: 



The fields e and h are introduced from the corresponding field strengths, i.e. Yiq = —a^Ci 
and Yij = —a?eijk ■ The gauge action is canonical in terms of E and B and not in terms of 
e and h. The evolution equations of the canonical modes derived from the action and their 
explicit form is: 



E 

7x 



J_d_ 



B 

(VXE), 



0, 



(A.2) 
(A.3) 
(A.4) 



where the possible presence of the Ohmic current has been included for completeness even 
if conducting initial conditions will not be considered explicitly. The system of Eqs. (IA.2p - 
(IA.4p . in the absence of electromagnetic sources, is invariant under the generalized duality 
transformation E — > —B, B ^ E and \/X — t- [371 |38]. The conventions for the Fourier 

transform are: 



Bi{x,T) 



(27r)3/2 



£kBi{k,T) e 



-ik-x 



Eiix, t) 



d^kEi{k,T)e-'^-\ (A.5) 



(27r)3/2 

Consequently the fluctuations of the magnetic and electric energy densities is given by: 

1 



(27r)3/2 Sna^ 
1 



J d'k ]^B,{k,T)B,{q- k,T) -^P^{k,T)S^^\q) 
J d'k T)E,iq- k,T)-^ PEik, t)5^'^ (q) 



The electric and magnetic anisotropic stresses are defined as 



. (A.6) 



nlf(g,r) 



nSf(g,r) 



d^k 



Aira'^ J (27r)3/2 
d^k 



J (27r)3/2 



B,{k, T)B,{q- k, t) - -^B^{k, T)B,^{q- k, r) 
E,ik, T)Ej{q- k, r) - -fE^{k, T)E^{q- k, r) 



• (A.7) 



The stochastic averages of the fluctuations variables defined in Eqs. (]A.6|) - (1A.7|) are all 
vanishing, i.e. using Eqs. fl2.5p -( l231) . {6pb{x,t)) = and {6pe{x,t)) = and similarly for 
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the anisotropic stresses. The second order correlations of the magnetic energy density and 
of the anisotropic stress are 

27r2 



(5pB(g,r)5pB(p,r)) = =^QB(g,r) 5^'\q + p), 

q6 



where 



QB(g,r) 



QBn(g, 



128 7r3 aS 

^3 



PBik,T)Pj,i\q-k\,r) 
\q-k\'^ 



^p{k,q), 



T 



288 773 a^ir) 



(Pk 



PBik,T)PBi\q-k\,T) 



k^ 



\q — k\^ 



An(fc,g). 



The functions Ap{k, q) and An(/c, q) are defined as 



Ap(fc,g) 
An(A;,g) 



1 + 



[fc • (g - fc)]^ 

[k-{q-k)f 
k'^\q-k\^ ' q 



(A.8) 
(A.9) 

(A.IO) 
(A.ll) 

(A.12) 



6 

+ - 



k ■ {q — k) 



[k ■ (g - 



+ 



9 

7 



k''\q-k\^ -2[k-{q-k)\f + 



k'^\q — fcp 
[k-{q-k)Y- 



(A.13) 



k'^\q — fcp 

The functions Ap{k,q) and An(A;,g) coincide for magnetic and electric degrees of freedom 
since both E and B are solenoidal fields: B is solenoidal because of the absence of magnetic 
monopoles while E is solenoidal because the pprotoinflationary plasma is globally neutral 
and any electric charge asymmetry is absent. The explicit expressions of the power spectra 
of Eqs. flA.lOp and flA.lip is obtained by using the power spectra of Eqs. (I2.2l) - fl2.3p . 

For typical wavelengths larger than the Hubble radius the second-order spectra including 
the slow roll corrections are given by: 



QB(g,r) = OB(g,e,/) 



a \9B{eJ) 



QBu{k,T) = OBu{q,e,f) 



a \9Bie,f) 



(A.14) 



where gsi^, f) = 4/ — 8 + 4e /. The fc-dependent amplitudes appearing in Eqs. (]A.14p are: 



OB{q,e,f) = H'CM,^)^B{f,e,q)(^ 

\qp 

0B(g, e, /) = CMf, e) £Bn(/, e, q) 



mBn(<:,/)-l 



(A.15) 



where mB(e, /) = rriBni^, /) = 9 — 4/(1 + e). The functions Cb(/, e) and CBuif, e) are given, 
respectively, by: 

94/(1+^) 4 

^^^r^[/(l + e) + 1/2], CBn(/,e) = -CB(/,e). 



CB(/,e) 



(A.16) 
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The functions >Cb(/, e, q) and jO-Buif, q) are: 

8[/(l + e) + l] 



'^B(/,e,g) 



2/(l+e)-4 



3[4/(l + e) - 5][4 - 2/(1 + e)] 3[4 - 2/(1 + e)] Vgo/ 



+ 



4/(l+€)-5 



'^^Bn(/, e, ?) 



5 - 4/(1 + e) Vgme 

2[17-2/(l + 6)] 

15[4/(1 + e) - 5][4 - 2/(1 + e)] 3[4 - 2/(1 + e)] \qoJ 

4/(l+e)-5 



(A.17) 



/ q ^ 2/(l+e)-4 



The comoving scale Qp = 0.002 Mpc~^ is the usual pivot scale at which the power spectra 
of the scalar curvature are assigned. The value of qo has been chosen 0.001 while ^max is 
related to the maximal amplified frequency of the magnetic field spectrum. 



B Evolution equations in the a parametrization. 

In the a-parametrization the Hamiltonian and the momentum constraints read, respectively, 



dh 
da 



2 2 

K a 



2(q; + 1)^ a 
dC 4 



MO 



MO 



da 



RJ^ + R^[l + i?b(a)]^7b + 7 



3 Vt 



cO 



aOr 



(B.l) 
(B.2) 



where R\){a) denote the baryon-to-photon ratio Rb{a) — 3(Qbo/^Mo) ck/(4-R7) — 0.215 a. 
The remaining two equations stemming from the perturbed Einstein equations can be written 
as: 

K^i 3 



d'^h 5a + 4 dh 
+ 



da'^ 2a{a + 1) da 2{a + 1] 



d'^Q ^ 5a + 4: dQ 



a2(« + 1) 
12 



R^6^ + Ri/6i/ + Ayf^B 



H — —rriRuf^v + R^cfb] 

a'^ia+l) 



da^ 2a{a+l)da 2(a + 1) a^{a+l) 
where Q— (/i + 6^); the evolution equations of the neutrinos obey. 



2dh 



3x/aTT 3aa' 



d9^ 
da 



d6^ ^ 
da 

da^ _ 

'da ~ IL^a + t 15 20 ^^TT' 

dJ^„t K 



K 



-5. 
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2 dQ 3 kT, 



u3 



da 2{2i+l)y/^n~T 
The evolution equations of the dark-matter sector obey instead 

d6^ _ idh de^ 

da 



2 VoTT 2 da 



da 



a 



(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 

(B.8) 
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In in the tight-coupling approximation the equations of the baryon-photon system are: 

d9^h RhO-yh _ k'^S^ At^ (Qb - 4(7b) g^ 

da a{R^ + 1) ~ 8A/^TT(i?b + 1) 8v/^r+T(i?b + 1) ' 
d5^^_2_d^ 2dh dS^ _ g-,b 1 dh 
da 3VaTT Sda' da 2^^TT 2 da' ^ ' ^ 
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